Superconducting phase qubit based on the Josephson oscillator with strong 

anharmonicity 
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We propose a superconducting phase qubit on the basis of the radio-frequency SQUID with the 
screening parameter value Pl = {2n /$o)LIc ~ 1, biased by a half flux quantum <l>e = 'l>o/2. 
Significant anharmonicity (> 30%) can be achieved in this system due to the interplay of the cosine 
Josephson potential and the parabolic magnetic-energy potential that ultimately leads to the quartic 
polynomial shape of the well. The two lowest eigenstates in this global minimum perfectly suit for 
the qubit which is insensitive to the charge variable, biased in the optimal point and allows an 
efficient dispersive readout. Moreover, the transition frequency in this qubit can be tuned within 
an appreciable range allowing variable qubit-qubit coupling. 

PACS numbers; 85.25.Dq, 74.50.-Fr, 03.67.Lx, 03.67.Pp 



The superconducting qubits based on the Josephson 
tunnel junctions (see, e.g., the reviews in Refs. [J 
have already demonstrated their great potential for the 
quantum computation [3]. The so-called phase qubits 
present the class of devices which are particularly suitable 
for integration with microwave on-chip transmission lines 
and resonators, i.e. the elements which significantly ex- 
tend the scope of the quantum circuit designs ^ . These 
qubits are based on the energy quantization in the shal- 
low wells of the inclined cosine Josephson potential Q. 
This shape is ensured either by finite bias current Is with 
the value slightly below the critical current of the Joseph- 
son junction Ic or a finite flux bias $e applied to the qubit 
loop (in the case of a loop configuration of the circuit) 
In both cases, the energy potential can be approximated 
by the cubic parabola with a smooth energy barrier iso- 
lating the well from one side and allowing escape out of 
this well enabling a simple readout. 

The low depth of the cubic parabola well leads to an- 
harmonicity, viz. successive reduction of the transition 
energies AEn = {En+i — En), n = 0, 1, from bottom 
to top, necessary for the qubit operation within the basis 
states |n = 0) and |n = 1), excluding unwanted excitation 
of the higher energy states {n > 1). Usually the phase 
qubit is designed such that for appropriate phase bias the 
cubic potential well includes three-four energy levels with 
anharmonicity of a few per cent 0, Q . This is achieved by 
adjusting the plasma frequency of the Josephson junction 
both by designing appropriate parameters of the junction 
and, possibly, by applying external capacitor shunting. 
The lowering of the energy barrier by applying the so- 
called measuring pulse, makes possible the reduction of 
the number of the levels to two (n = and 1), with no- 
tably different rates of escape to a running-phase state (in 
the case of current bias), or to the lower-energy state in 
the adjacent well (in the case of the loop configuration of 
the qubit). The large (but finite) difference of these tun- 
neling rates sets the maximum theoretical value for the 
fidelity of such measurement to 96.6%. In the carefully 



designed and optimally biased qubit the best experimen- 
tal fidelity values approach 90% . The main disadvan- 
tage of such phase qubit is the necessity of resetting it 
after each measurement. 

In contrast to the charge fsl , charge-phase [l^ , flux 
nj, transmon p^] and recently proposed, the so-called 
fluxonium [ist qubits, the conventional phase qubits can- 
not inherently operate in an optimal point, i.e. in the 
symmetric working point insensitive in the first order 
to the noise that could give drastic improvement to the 
qubit performance [^, (The exceptions are the re- 

cently proposed three-junction interferometer circuit 15 1 
and the so-called camelback potential phase qubit based 
on the two-junction SQUID [161.) Moreover, the limited 
anharmonicity of the phase qubit makes the observable 
reactance impedance (i.e. the Josephson inductance) val- 
ues in the ground and exited states hardly distinguish- 
able. This poses serious problems for dispersive read- 
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FIG. 1: (a) Electric diagram of the qubit coupled to a resonant 
circuit and (b) possible equivalent compound (two-junction 
SQUID) circuit of the Josephson element included into the 
qubit loop. Capacitance C includes both the self-capacitance 
of the junction and the external capacitance. Due to inclu- 
sion in the resonant circuit of a Josephson junction JJ', the 
resonator may operate in the nonlinear regime, enabling a 
bifurcation-based readout. 
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out schemes, which proved advantageous where apphca- 
ble [13, U, 121 and have allowed quantum nondemoli- 
tion measurements as well as hi gh fi delity measurements 
based on bifurcation amplifiers [iSl . fioj . 

In this paper we propose an improved phase qubit 
with significant anharmonicity, in which the manipula- 
tion and dispersive readout are both possible in a sym- 
metry point. The circuit diagram of our qubit is shown 
in Fig. la. It comprises the superconducting loop with 
geometrical inductance L closed by a Josephson junc- 
tion, generally shunted by an external capacitance, and 
rf-driven LuCrGr resonance circuit inductively cou- 
pled to the qubit loop. The peculiarity of this qubit 
is the unity value of the SQUID screening parameter 
/3i = (27r/$o)i^c ~ 1, where $o = h/2e is the flux quan- 
tum. This can be achieved by an accurate design of the 
circuit including replacement of the single junction with 
a two-junction SQUID, allowing more precise adjustment 
of the resulting critical current Ic (see Fig. lb). 

The potential energy of the stand-alone qubit biased 
by external magnetic flux <I>e includes the magnetic and 
Josephson components and can be written as 

[/((/), 4>e)^EL [0.5((/. - (t>ef - /3l(1 + COS (/.)] , (1] 

where El — {^o/2tt)'^/L = Ej/Pl is characteristic mag- 
netic energy associated with the loop inductance, the 
Josephson couphng energy Ej = (<I>o/27r)/c, the phase 
variable (p and the phase bias (pe — 27r$e/'&o- For small 
values oi (3l «C 1, the potential Eq. ([T]) yields the almost 
parabolic shape of the global single well (first term in 
Eq. Il])), whereas for the values Pl appreciably greatei 
that 1, the series of wells are superimposed on the globa' 
parabola, so the bottom parts of these local minima car 
also be approximated by the quadratic parabolas. In the 
case of large density of the levels within these parabolas 
the energy spectrum is also close to that of a harmonic 
oscillator. So, neither of these cases allows significant an- 
harmonicity necessary for the efficient qubit operation. 

The essentially different shape of the potential Eq. ^ 
with — t: is, however, achieved for /3i w 1, i.e. when 
the quadratic magnetic term is partially compensated by 
the quadratic term in the Josephson energy expansion 
near the bottom {U — 0) of the single well centered at 
(y9 = — TT = 0, i.e. 



U{ip) « Ej 



(Pl - 1) 
2/3l 



1 



(2) 



In the ultimate case oi Pl = 1, Eq- © provides the ob- 
tuse shape of the quartic parabola. Taking into account 
the finite kinetic energy of the system in correspond- 
ing Schrodinger equation, Q^/2C = —4Ecd^^, where 
Q = —i{2e)d(j, is the charge operator [i^, one can per- 
form quantization of the system. Application of the qua- 
siclassical quantization rule of Bohr-Sommerfeld yields 
for this quartic oscillator the energy levels obeying the 



4/3 power law [2J|: 



(3) 



with prcfactor e which in terms of the parameters of our 
circuit is equal to 

e = 2-5/33 [^/K(l/2)]4/3(£;^i?f )i/3 ^ i,Q{^EjElY/\ 

(4) 

where K(A:) is the complete elliptic integral of the first 
kind. Thus, the energy spectrum in the quartic potential 
takes intermediate position between the equidistant spec- 
trum of the harmonic oscillator En cx (n -I- 1/2) and that 
of the rectangular well, En oc (n -I- 1)^, having extremely 
high anharmonicity. Expressions ([3]) and (|4]) are exact for 
the higher levels (n ^ 1) and large "mass" (capacitance 
C) , ensuring the very large ratio of the Josephson energy 
Ej to the charging energy Ec — e^/2C. An estimate of 
the anharmonicity factor in this quasiclassical approxi- 
mation can be immediately obtained from Eq. Q : 



5qc = (Ai;i - AEo)/AEo « 26%. 



(5) 
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FIG. 2: (Color online) Position of the lowest six levels (solid 
lines) in the potential Eq. ([1} for 0e = tt as a function of 
parameter j3L for typical values of L and C, yielding Ej /Ec ~ 
El/Ec « 5.1 X 10^. With an increase of Pl, the spectrum 
crosses over from that of the harmonic oscillator type (left 
inset) to the set of the doublets (right inset), corresponding to 
the weak coupling of the oscillator-type states in two separate 
wells. The spectrum in the central region /3l « 1 is strongly 
anharmonic. The dashed line shows the bottom energy of the 
potential U{(j),(f)e = tt), which in the case of /3l > 1 is equal 
to ~AU ^ -l.hELiPh - if/PL (in other words, AU is the 
height of the energy barrier in the right inset) [20l. [2H. The 
dotted (zero-level) line indicates the energy in the symmetry 
point (f} — 0, i.e. at the bottom of the single well (/3l < 1) 
or at the top of the energy barrier (Pl > 1). The black dot 
shows the critical value PI at which the ground state energy 
level touches the top of the barrier separating the two wells. 
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The numerical solution of the corresponding Schrodinger 
equation with potential energy Eq. (1) yields in the limit 
Ej / Ec ^ 1 an even larger value of the anharmonicity fac- 
tor, 5 « 33% (see the energy spectrum in Fig. 2). These 
values substantially exceed the typical anharmonicity 
values of the conventional phase qubit, |i5phasc| ~ 3%, 
for the number of levels inside the cubic-parabola well 
equal to four 0, and transmon-qubit, l^transmoni_~ 
{EJ?,EjY/^ < 5% for optimum values Ej/E^ > 50 
Moreover, in contrast to the negative values of S in these 
examples, the series of the energy levels in the quartic po- 
tential has positive S > 0, i.e. corresponds to successively 
increasing level spacings AEi < AE2 < AE^... 

Such a large, positive anharmonicity is a great advan- 
tage of the quartic potential qubit allowing manipula- 
tion within the two basis qubit states |0) and |1) not 
only when applying resonant microwave field, i^^w ~ I'm, 
but also when applying control microwave signals with 
large frequency detuning or using rather wide-spectrum 
rectangular-pulse control signals. The characteristic 
qubit frequency via = AEo/h and the anharmonicity 
factor S computed from the Schrodinger equation for the 
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FIG. 3: (Color online) (a) The qubit frequency as a func- 
tion of parameter /3l for fixed L = 50 pH and several val- 
ues of capacitance C — 0.1, 0.3, 1.0 and 3.0 pF (from top to 
bottom), corresponding to the values of the ratio El/Ec ~ 
1.7 X 10*, 5.1 X 10**, 1.7 X 10^ and 5.1 x 10^ (b) Anharmonic- 
ity parameter 5 as a function of parameter I3l for the same 
as in (a) inductance L and capacitance values (from top to 
bottom). 



original potential Eq. ([T]) in the range 0.9 < /3l < 1-02 are 
shown in Fig.[3l One can see that the significant range in 
the tuning of the qubit frequency within the range of suf- 
ficiently large anharmonicity (^ 20 — 50%) is attained at 
a rather fine (typically ±1-2%) tuning of /3l around the 
value Pl — ^- Such tuning of /3l is possible in the circuit 
having the compound configuration shown in Fig. lb. For 
values of /3l > 1, the symmetric energy potential has two 
minima and a barrier between them. The position of the 
ground state level depends on Pl and the ratio of the 
characteristic energies Ej/Ec = PlEl/Ec- The value of 
/?L at which the ground state level touches the top of 
the barrier sets the upper limit /3£ for the quartic qubit 
(marked in Fig.[2]by sohd dot). At Pl > /3£, the qubit 
energy dramatically decreases and the qubit states are 
nearly the symmetric and antisymmetric combinations 
of the states inside the two wells (see the right inset in 
Fig. 2). Although the qubit with such parameters has 
very large anharmonicity and can be nicely controlled by 
dc flux pulses [U [i^l , its readout can hardly be accom- 
plished in a dispersive fashion. 

Another advantage of the phase qubit having the en- 
ergy potential of the shape close to the quartic one is a 
strong dependence of its Josephson inductance Lj($e, n) 
on the quantum state |n). The observed value of the re- 
verse inductance is related to the local curvature of the 
dependence of corresponding energy En on flux <i>e (see, 
e.g., Ref.[il), 



9$? 



(6) 



where / is the operator of supercurrent circulating in the 
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FIG. 4: (Color online) The values of the Josephson induc- 
tance of the quartic potential qubit in the ground (solid lines) 
and excited (dashed lines) states calculated for the geometric 
inductance value L — 50 pH and the set of capacitances C, 
increasing from top to bottom for both groups of curves. 
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qubit loop. The dependence of the reverse inductance 
Lj{^e = $0/2, n) calculated numerically in the two low- 
est quantum states (n = and 1) for L = 50 pH and the 
same set of capacitances C as in Fig. [3] is shown in Fig.lD 
One can see that the ratio of the geometrical to Joseph- 
son inductances L/Lj takes large and very different val- 
ues that can be favorably used for the dispersive readout, 
ensuring a sufficiently large output signal. Note that for 
/3l < 1, both inductances Lj{n — 0) and Lj{n — 1) are 
negative, whereas at /3l > 1 the inductance Lj{n = 1) 
changes the sign to positive. 

The readout of this qubit is based on the measurements 
of the reactive part (inductance) of the loop impedance 
probed by a low- frequency ac signal, / <C lAnoi suffi- 
ciently small amplitude [25( (see also Ref. [lO|). This 
signal is supplied by a rf-driven oscillator (Fig. la) as 
an alternating biasing flux, <I>e = 0.5$o + MIl or (j)^ — 
TT + 6(j)e, where S4>e — acos(27r/i) with a ^ 1. Here 
M = K{LLfiY^'^ is mutual inductance, k, a dimcnsionless 
coupling coefficient and II, the ac current in inductance 
Lji. Coupling of the qubit to the resonance tank circuit 
causes renormalization of the circuit inductance (see, e.g. 

EEl), 



= Ln{l - K^L/Lj{n)), 



(7) 



and the resonance frequencies Un = [L^^C]~^/^, where 
n = and 1. The relative difference of the resonance 
frequencies for the qubit in the excited and ground states 
is 



Slu 



LOi - ^0 



' l^K^L/Lj(0) 
1 - K^L/Lj{l) 



1. 



(8) 



Figure [H] shows this relative frequency shift versus pa- 
rameter (3l ■ One can see that for the rather conservative 
value of dimcnsionless coupling k = 0.05, the relative 
frequency shift can achieve the easily measured values 
of about 10%. The efficiency of the dispersive readout 
can be improved in the non-linear regime with bifurca- 
tion 18| . With our device this regime can be achieved in 
the resonance circuit including, for example, a Josephson 
junction (marked in the diagram in Fig. 1 by a dashed 
cross). Due to the high sensitivity of the amplitude 
(phase) bifurcation to the threshold determined by the ef- 
fective resonance frequency of the circuit, one can expect 
a readout with high fidelity even at a rather weak cou- 
pling of the qubit and the resonator (compare with the 
readout of quantronium in Ref. (26|). Further improve- 
ment of the readout can be achieved in the QED-based 
circuit including this qubit [27| . 

The loop configuration and frequency detuning of the 
quartic qubits should allow their inductive coupling with 
variable strength keeping both qubits in optimal points. 
Variable coupling of the optimally biased qubit to a su- 
perconducting resonator is also possible. More sophis- 
ticated coupling of the pairs of quartic qubits can be 



0.1 



ffl 



0.01 
0.97 



L = 50 pH 


C = 3.0pF / 

C = 1.0pF / 

C = 0.3pF / , 

c = o.ipF y'^ 







0.98 



0.99 
Pl 



1.00 



1.01 



FIG. 5; (Color online) The resonance frequency shift in the 
circuit due to excitation of the qubit with the inductance value 
1/ = 50 pH and the set of capacitances C, decreasing from top 
to bottom. The dimensionless coupling coefficient k = 0.05. 



accomplished, for example, using a Josephson-junction 
coupler in a fashion recently proposed by Harris et al. 

In conclusion, we have shown that the phase qubit of 
the rf-SQUID configuration with parameter /3l ~ 1 and 
flux bias 'I'e = <I'o/2 has remarkable characteristics. Still, 
we expect that implementation of this qubit requires the 
solution of several experimental problems. For exam- 
ple, due to a high sensitivity of the qubit parameters to 
the magnitude of /3i, whose optimum values lie within 
a rather narrow range (±1-2%), particular precaution 
should be taken against fluctuations in the line control- 
ling the effective Josephson coupling in the circuit (see 
Fig. lb), because otherwise it may cause significant de- 
phasing of the qubit. Furthermore, flux bias <I>e — 'I'o/S 
should also be set as precisely as possible. Experimen- 
tally, it can be realized either by freezing the flux $e = 'i'o 
in the main loop having a symmetric gradiometer config- 
uration (29I I , or by including in the loop of a Josephson 
TT-junction '305 with a sufficiently high critical current 
ensuring the steady phase shift of tt. 

Of course, similar to properties of the conventional 
types of the phase qubits, the coherence characteristics of 
the quartic qubit will be strongly dependent on the ma- 
terial properties of the circuit. Minimizing the losses due 
to the qubit coupling to microscopic degrees of freedom 
(two-level systems) inside the dielectrics surrounding the 
superconducting circuit (the substrate, insulator inside 
the capacitor, the junction barriers, etc.) play crucial 
role for improving the qubit coherence [3| • Since the op- 
eration and tuning of the quartic qubit is possible with- 
out leaving the optimal point, one may expect a weaker 
coupling of the qubit to these microscopic two-level sys- 
tems located inside dielectrics and, therefore, a better 
quantum coherence. Moreover, the zero persistent su- 
percurrent circulating in the qubit loop at the optimum 



5 



bias, (pe = TT, may also reduce the effect of quasiparticle 
tunneling on the qubit coherence. Probably, such weak- 
ening of the qubit coupling to external degrees of free- 
dom can explain reasonably good coherence characteris- 
tics (TRabi ~ 60 ns) of the Nb camelback qubit operated 
in an optimal point at zero persistent current [l6|. Any- 
way, the properties of the proposed quartic qubit will be 
clarified in experiment which is currently in the prepara- 
tion stage. 
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